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Abstract — In the Markov and covariance interpolation prob- 
lem a transfer function W is sought that match the first 
coefficients in the expansion of W around zero and the first 
coefficients of the Laurent expansion of the corresponding 
spectral density WW*. Here we solve an interpolation problem 
where the matched parameters are the coefficients of expansions 
of W and WW* around various points in the disc. The solution 
is derived using input-to-state filters and is determined by 
simple calculations such as solving Lyapunov equations and 
generalized eigenvalue problems. 

I. INTRODUCTION 

The problem of designing filters from covariances and 
Markov parameters has been studied before in numerous 
papers [17], [13], [14], [15], [16], [20], [21], [18]. Skel- 
ton et. al. call a stable model matching Markov param- 
eters Hq, Hi, ... , Hq^i and covariances Rq, Ri, . . . , Rq-i 
a q-Markov COVariance Equivalent Realization (q-Markov 
COVER) and they have shown that if the data satisfies a 
particular consistency condition (which can be avoided using 
a variable input variance as in [4]), there are many such 
q-Markov COVERs and they are parameterized by a set 
of unitary matrices. One of the parameters considered "as 
known" in the classical q-Markov COVER theory is the 
variance of the input noise. In [4], [5] the author proposed 
a method for designing minimal degree realizations using 
the variance of the input noise as a design parameter which 
enabled a realization of lower degree to be determined. In 
fact, that method guarantees a unique stable solution for 
generic data. Here, using input-to-state filters, we solve an 
interpolation problem where the matched parameters are the 
coefficients of expansions of W and WW* around various 
points in the disc. We could for example consider matching 
the constraints 



noise 



91, 



for some points pi , • • • , p„ in the unit disc and similarly for 
WW* . Most results in [4] carry over to this more general 
problem. A formal definition of the problem considered is 
given in the next section. Another approach to this problem 
was taken in [19]. The main objective of that paper was to 
prove existence of a fixed point for the Steiglitz-McBride 
algorithm and a different kind of normalization was used. 

II. Problem formulation 

First the Markov and Covariance interpolation problem as 
formulated in [4] is described and then input-to-state filters 
are introduced for treating the generalized problem. 
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Fig. 1. System considered 



A. The Markov and Covariance interpolation problem 

We consider a SISO system where a deterministic control 
signal V and a stochastic noise signal w are fed through the 
same system W to produce the output y as depicted in Fig.[T| 

Define u ^ v^w and let v be the control input and w an 
additive noise term. Assuming that the transfer function W 
is rational and of McMillan degree 77, it can be described by 
a minimal state space system 



Xj+i = Axj 



-Bu 
Vu 



(1) 



where A £ C'^", B e C"^\ C £ C^^" and V e C. The 
output y is the superposition of the outputs due to each of 
the inputs v and w. Therefore, data from the system can be 
obtained by the following idealized experiments - or in any 
other practically more suitable way. 

First, determine the output when the noise w is zero and v 
is a unit impulse, yielding the Markov parameters (impulse 
response parameters) 



Ho, Hi, 



,Hf 



(2) 



Second, determine the output when the control v is zero 
and w is mean zero white noise with unknown variance A G 
C. Assuming W is asymptotically stable then this system 
provides a realization of a stationary stochastic process, and 
by truncated ergodic sums the covariances 



Rq, Ri, . . . ,Ri 
can be estimated such that the condition 
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(4) 



is satisfied. 



B. Input-to-State filters and interpolation 

In order to analyze a signal it is useful to consider a new 
signal obtained by applying an input-to-state filter [10], i.e. 
if Uk is our original signal we define the new state vector 
by 

Xk = Axk^i + Byk, X(i = 0, (5) 

where e C"''", B e C"""^ and the eigenvalues of A lies 
in the open unit disc. Note that the state x in O is not the 
same as the new "artificial" state x defined in Q from the 
"inputs" j/fe. 

Consider the input to state map G: 

G{z) ■={1 - zA)-^B (6) 
where we will assume that (A, B) is a reachable pair, i.e. 

T ^[ B AB ... A'^'^B ] (7) 

is full rank. 

A wide class of interpolation problems can now be ap- 
proached in a unified framework by expressing the interpo- 
lation constraints as inner products with the input-to-state 
map G. Let (•, •) denote the standard L2 inner product on 
the circle, and for vector- and matrix-valued functions Fi 
and F2 define 

1 



{Fi,F2) = 



2ti 



Fi{e'")F*{e"')d0, 



1 



where the integral is evaluated elementwise and F* denotes 
the adjoint of G, i.e. 

F*{z)^F{-z-^r, 

where the superscript * denotes the usual complex conjugate. 
Notice that we will allow the L2 inner product between 
two matrix-valued functions Fi and F2, possibly of different 
sizes, provided that the product Fi{z)F2 (z) is well defined. 

In the special case where A — diag(pi,p2, • ■ • ,Pn), \Pj\ < 
1, and B ~ [1, 1, • • • , 1]^, the scalar function on the fc:th row 
of G(z) is 

9k{z) = -r 

1 ~PiZ 

Then from the Cauchy's integral formula {f,gk) = fip*), 
i.e. the values of / at the selected points can be expressed 
in terms of the inner product. 

In the other special case where all interpolation points are 
at the origin, i.e. the values of the function and its derivatives 
at zero are interpolated as in the Caratheodory interpolation 
problem, then we could chose 



(8) 
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so that 



G(z) = 



(9) 



The states are then the n most recent outputs and it is easy 
to see that the covariance of the state is a Toeplitz matrix as 
in (IHi. 

In practice one could be interested in having a mixture of 
interpolation conditions on the function values at different 
points and on some of its derivatives, and this can be 
accomplished by considering for example ^-matrices with 
some particular Jordan structure. To be able to find a B 
such that {A, B) is reachable it is necessary that A is cyclic, 
so there can not be more than one Jordan block for each 
interpolation point (eigenvalue of A). 

Now given some G, if is a matricial spectral measure 
of the input (i.e. y the input to G) the state covariance S 
will satisfy [10] 



G{e^^)dfi{d)G{e 



(10) 



For the more general input-to-state filter it is more difficult 
to know what is the structure of the state-covariance matrix. 
In Theorem |2] below, a result from [10] describing the 
feasible structures is stated, but first we need to remind the 
reader of a well-known result. 

Lemma 1: The matrix Q defined by 



A 



Q = {G,G), 



(11) 



is the Reachability Gramian solving the discrete time Lya- 
punov equation 



g ^ AQA* + BB*. 



(12) 



Since {A, B) is assumed to be a reachable pair, Q is 
invertible. 

Proof: Note first that 



AG{z) = z-^{G{z)~B), 



(13) 



and then multiply ( fTTT i with A from the left and A* from the 
right to obtain 

AgA* = {AGG*A*,l) 

= {z~\G - B)z{G* - B*),l) 

= {GG* + BB* - BG* - GB* , 1) 

= g + BB* - BB* - BB*. (14) 

The last step follows by observing that G is analytic in 
the unit disc and thus (G, 1) = G(0) = B, and similarly 
(G*,l) = G*{oo) = B*. 

Since A is assumed to be asymptotically stable the solution 
to the Lyapunov equation is unique, and this completes the 
proof. ■ 

Theorem 2: A positive definite matrix E is a state- 
covariance matrix for a suitable input process if and only 
if it is of the form 

j:^^{Mg + gM*) 

for a matrix which commutes with A. Furthermore, 
any such matrix A4 is uniquely defined modulo an additive 
imaginary constant al with a G jM. 



Another way to describe the structure of the state covari- 
ance E is that it satisfies the equation [8] 

S - AT,A* = BL + L*B* 

for some L. 

Let 'H2 denote the Hardy space of functions that are 
analytic in the unit disc with square-integrable radial limits, 
and define 



K^'H2eb{z)H2, 



(15) 



where b[z) = det(z/ - A*)/det(/ - zA) is a Blaschke 
product with poles at the eigenvalues of A. In fact, h{z) is 
the inner, or Douglas-Shapiro-Shields, factor of G{z). Then 
JC contains all functions in T-L2 which are orthogonal to those 
that vanish on the spectrum of A*, and it is usually called 
the coinvariant subspace. By [9, Prop. 4] the elements of 
G{z) form a basis for JC, so any / e /C can be written 
/(z) — CG{z) for some vector C, and then 

det{I-z{A-BC))-det{I-zA) ^ ^ 
^(") = det(/ - zA) ^'^^ 

We also need to take inner products between elements in 
K. and ■H2, and then the following formulas are useful. 
Lemma 3: If f{z) G H2 then 

(/,G) = B*f{A*) (17) 

and 



(GJ) = f{A)B. 



(18) 



Furthermore, 



(G,/G) = f{A)g. (19) 
Note that it is important here that / is a scalar function. 
Proof: Since / e -H2 and G e H^'"'^ they have series 
expansions 



and 



Then 



00 00 

G(z) ^^Gfcz'^ = ^^'=Sz'= 



fe=0 



fe=0 



00 00 



k=Q e=o k=0 

and the formula for (/, G) follows by considering the com- 
plex conjugate. 
Finally, 



00 00 00 



(G,/G) = (A'=i3z^/,z^G™z™) 

00 00 



£^0 m^O 

00 00 

Em^ E ^'"^^*(^" 



which concludes the proof. ■ 
Remark 4: Alternatively, this could be proven by consid- 
ering generalized Cauchy kernels 



as in [10]. 

Estimation of the parameters from data can be performed 
by applying the input-to-state filter and then using standard 
techniques, see [2] for examples of filter bank data analysis. 

III. The global optimization problem 

We will assume here that the spectral measure in (fTot is 
given by 



where 



(20) 



fe=0 



i.e. is analytic in the unit disc (so the sum converges for all 
z in the unit disc), and in this class of spectral measures 
we will find the one allowing the maximal input variance A 
meanwhile satisfying the following interpolation conditions: 



{GWKW*G\1) = S, 



(21) 



where the state covariance S satisfies the condition in 
Theorem 121 and 

(G, W) = H, (22) 

for an arbitrary nonzero state-Markov vector H. 

The interpolation constraint in (I2II was considered in, 
for example, [8]. The interpolation constraint in (l22t can be 
recognized as a special case of the Lagrange-Sylvester inter- 
polation as studied in [1, section 16]. Here, both constraints 
are enforced simultaneously. 

Thus the optimization problem considered is: 



(6) 



max 
W eH2 

A e IR+ 

s.t. 



A, 



(GWAW*G*,1) 
(G, W) = H. 



Let S be an (n x n) Hermitian matrix and C be an (1 x n) 
vector consisting of Lagrange multipliers, the Lagrangian is 
then 

C{W,A) = A + tr{(E- (GWAW*G*,1))2} 
+a{G,W)-H). 

We can rewrite it in the following form 

C{W,A) = A + tr {ES} - (T^AT^^, G*SG) 
-(H + (CG, W) . 

where 

G*EG ^B*{I- z-'^A*y^E{I - zA)-^B 



and 



QG^C{I-zA)-^B 



are scalar functions. 

Before taking the maximum we write it in the form 

C{W,K) = (CG, W^) + ((1 - W^*G*SGiy)A, 1} 
+tr{SS} - CH 

Note: Sup{£(PF,A)|PF G H2,A > 0} < oo only if G*SG 
is in the "positive cone", i.e. it is non-negative for all z on 
the unit circle, and 

{W*G*^GW,l) >l. (23) 
Maximizing over A while assuming ( l23T l it must hold that 
A (1 - Ty*G*SGH/, 1) = 0, (24) 
and since A 7^ 0, equality must hold in (l23T l. i.e. 

(T^*G*SGM/,1) = 1. (25) 

Maximizing over W shows that the following variation 
has to be zero for all 5W 



Therefore 



(CG - 2AG*5GVF, 5W) = 0. 



CG = 2AG*SGM^ + V* 



(26) 



(27) 



where F £ and V{Q) = 0. From ^ the poles of V* 
has to be poles of G*. Furthermore, V ^ JC follows by 
considering the partial fraction expansion of (l27t . so there 
must be a vector v such that 



V{z) = vG{z) 
Then the transfer function W will be given by 



(28) 



W ■ 



1 



(G*SG)-i(CG-F*)A-i 



1 



= —{G*^G)-\CG-G*v*). 

Lemma 5: If S is non-negative we can factor G*'E.G as 

G*SG=(CGr(CG), (29) 

where ^ is a row-vector. 

Proof: Since S is non-negative and Hermitian it can be 
factorized as 



s = [ C2 



c ] 



en 



Then 



G*SG = ^G*efcaG 



(30) 



k=l 



and G*'E.G is a sum of elements in /C U /C*, where k* G /C* 
if fc G /C. Since all the terms in (|30] | are positive, by spectral 
factorization a vector ^ such that the sum is equal to G*^*£,G 
can be found. ■ 



For W to be analytic outside the unit disc it is necessary 
that the factor (^G)^* is cancelled, i.e. we need that 

2A(eG)(CG)'^ CG- 
From (|27li and it follows that 

CG = 2AG*<e VG + G*i^* . (3 1) 

Using (|25T l and dSTT l the dual function is 

^{^,<j) = ^cG,|^^+tr{i]ra-ci^ 

= 2Acr^;(T* + - C-H^ 

since (V^*, W) = and where Q was defined in ( fTTT l. 

To determine the last term (^i/, multiply OTl i with G* and 
integrate to obtain: 

C (GG^ 1) = 2A {G*CcrGG*, 1) + (y*G*, 1) 

the last term is zero and then 

(H = 2Aa{G{G*g-^H)G*,l)C = ^AaHC, 

where 

n = {G{G*g-^H)G\l) . (32) 

Lemma 6: The matrix T-L defined by ( l32b is the unique 
solution to the Stein equation 

n = AHA* + HB*. (33) 
Proof: As in the proof of Lemma [T] note that (13[ holds 
and then multiply ( l32l ) with A from the left and A* from the 
right to obtain 

AHA* = {AGG*A*{G*g-^H),l) 

= {z-'^iG - B)z{G* - B*){G*g-'^H),l) 

= ((GG* + BB* - BG* - GB*){G*g-^H), l) 

^n + B {{B* - G*){G*g-^H), 1) 

^{GB*{G*g-^H),l) (34) 

The second term in (|34| | is zero since the integrand is analytic 
outside the unit circle and G*(oo) — B*. 

The third term in (l34t is HB*, which follows by consid- 
ering the action on an arbitrary vector v; 

{GB*{G*g-^H),l)v = {GB*v{G*g-^H),l) 
= {GG*g-^H,l)B*v 
= HB*v. 

Since A is assumed to be asymptotically stable the solution 
to the Stein equation is unique, and this completes the proof. 

■ 

Remark 7: Note that if 

h{z) := [I - zA)-^H (35) 

then (/i, G) solves the Stein equation ( l33T l. and then by 
uniqueness (compare [9, Eq. (40)]) 

-H = (/i,G). 



Now, the dual optimality function is 



2Ag 





'2AH 





■ cr* " 







(36) 



Maximizing this expression over positive A 

ifii^, <j) = 2A{crgcr* - aUC) + (^C 

it is clear that {aQa* —aH(,*) has to be zero. In fact, if it is 
negative the optimal value of A would be zero and we have 
assumed that it is positive, and if it is positive the optimal 
value of A would be infinite. Furthermore, the following 
holds: 

Lemma 8: Given that W{z) = £ 'H2, the constraint 

(G, W) =H implies that TiC = Qa* . 



Proof: We know that (G, W) = H, and thus 



is a scalar function, so 

GG*c (^G*g-^ (g, ^ 



{GG*,^G)g-UG, 



aG_ 

Jg 

(G,(eG)G)g-i^G,^. 

m{A))gg-^(^G{A))-^ga* 
g<T* 



G)a* 



where we have used (fT9] l twice. 
The complementarity condition 



can be formulated as 



(jga* 



and then the dual problem is 

(D) 



min ^S^* 



s.t. 



ga* = HC 
aga* ^ 1 



where a and ^ are related by the Markov interpolation 
conditions. 

The variable A was eliminated above, but it can recovered 
by considering the dual of the dual. Let A be the Lagrange 
multiplier and use g<j* — to eliminate a 

L = - A((7g(T* - 1) 

= ^T^c - A^n*g^^nc + ^ 

= ^{E-A'H*g-^'H)C + ^ 

which leads us to maximize A as E — AUg^^H* is non- 
negative, i.e. 



max A, 
A 



s.t. T.- Awg-^n^Q 



The optimal A is now given by the largest positive value 
such that E — AH*g^^'H is non-negative definite, i.e. the 
smallest generalized eigenvalue of {Y.^'H*g^^T-L). 

Theorem 9: Given a state covariance E satisfying the 
condition in Theorem |2l and an arbitrary nonzero state- 
Markov vector H. Then, an optimizer W to problem (6) is 
given by W{z) — {(tG{z)) / {(^G{z)), where ^ is a nonzero 
solution to the equation 

{j:-AH*g-^n) e = o, 

A is the smallest generalized eigenvalue of {T,,'H*g~^'H), 
and finally a is determined by 

a* = g-^nC- 

Furthermore, if the smallest generalized eigenvalue has mul- 
tiplicity one the optimizer W is unique. 

Proof: This follows from the derivation above ■ 
As in the Markov and Covariance interpolation problem 
there is a special choice of the Markov parameters that reduce 
the problem to the equivalent of the Pisarenko method [11]. 
Namely, chosing H ^ B, then H — g and maximizing A 
under the constraint 

E - An*g^^n = e - Ag ^ o 

makes the corresponding W an inner function [11]. 

IV. Model reduction example 

To illustrate the method proposed here, a model reduction 
application is considered. The method proposed here is a 
generalization of the q-Markov COVER methods, that were 
initially proposed to be used for model reduction [14]. 

The transfer function from input 2 to output 1 of a 
portable CD-player is considered. This model, of order 120, 
is provided by SLICOT [3], and has been used by, for 
example, [12], [7]. The magnitude of the transfer function is 
depicted with a thick solid green line in Figure |2] There is 
a wide range of frequencies over which there are interesting 
features of the Bode plot. 

The given transfer function is a continuous time stable 
function. The bilinear map 

_ 1 - sT/2 



l + sr/2' 



(37) 



where T — 1/250, is used to transform the continuous time 
model into a discrete time model. 

Using three different input-to-state filters, reduced order 
models of degree 13 are designed. Our aim will not be 
to find the optimal interpolation point locations for this 
particular model, but to illustrate the way this choice effects 
the solutions. 

First an input-to-state filter as in (|8) was applied, corre- 
sponding to the Markov and covariance interpolation prob- 
lem described in section III-AI and the magnitude plot of 
the resulting model is depicted with a blue dashed line in 
Figure |2] 

Then, an input-to-state filter with 14 poles spread evenly 
around a circle with radius 0.95 was applied. The magnitude 
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Fig. 2. Model reduction 



plot, depicted in Figure |2] with a solid red line, is similar to 
the first one, but with a slightly smaller error for low and 
high frequencies. 

Finally, an input-to-state filter with 14 poles spread un- 
evenly around a circle with radius 0.9 was applied. The 
spread in frequencies were chosen to correpond to a log- 
arithmic spread in the frequency interval 10^ to 10^. In 
the discrete domain, the interpolation point locations are 
depicted with black plusses in Figure |3] together with the 
interpolation points of the two other filters. This choice of 
interpolation points is made to compensate for the frequency 
warping caused by the bilinear map (|37] |. The magnitude plot, 
depicted in Figure |2] with a black dashed-dotted line, shows 
an improvement of the fit in the frequency range where the 
poles were chosen. 

For comparison, a model of degree 13 is determined using 
a standard balanced truncation model reduction method and 
its magnitude plot is depicted in Figure |4] A good fit for 
the interval of frequencies where the magnitude is large is 
obtained. It is well known that weights can be applied to 
improve the fit for certain frequency regions. The choice of 
these weights, as well as the choice of interpolation points 
in our approach, should be made with the prior knowledge 
and requirements of the low order model in mind. 

V. Useful formulas for the user 

In this section we give simplified formulas for calculating 
the transfer function and state space representations of the 
W parameterized by ^ and (. It is also shown how to 
determine the state-covariances and state-Markov parameters 
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Fig. 3. Interpolation point locations 

from these representations of W. (Note that the problem 
considered in this paper is the inverse of determining the 
interpolation parameters from the model W.) These formulas 
will be important for applying the method proposed here. 

A. For transfer functions 

Since W = aG/£,G is a quotient between two functions 
in /C, it follows from ( fTSI l that it can be written as a quotient 
of two polynomials 

b{z) 



W{z) 



a{z) 



(38) 
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Fig. 4. Model reduction 



where 



biz) = det(/ - z{A ~ Ba*)) - det(J - zA) 



and 



a{z) = det(/ - z{A - EC)) - det(/ - zA). 



(39) 



(40) 



Clearly, z = is a zero of both a and h so it is cancelled 
out, which leaves a.W of degree at most n — 1. 

To check if the resulting transfer function W satisfies the 
interpolation conditions it is convenient to use (fTSl l to obtain 

(G, W) = W{A)B = b{A)a{A)-^B. 

To determine the state-covariance S corresponding to a 
particular W we can use the following formula from [10] 

S = K{GW,GW) = ]^K{^g + g^*) (41) 

where 5* = f{A) and / is the positive real part of WW* . 
The function / satisfies 

/ + /* = WW\ 

and it is clear that f — d/a, where a is given by ( |40] | and d 
solves the equation 



h{z)b*{z) = d{z)a*{z) + a{z)d*{z) 



(42) 



This equation has a unique solution d such that all the roots 
of d are outside the unit disc provided that all the roots of 
a also are outside the unit disc [22], [23]. 



Given a state space representation of W as in ([U the 
product GW has a realization 



(44) 
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Then the state-Markov parameter is given by 

(G, W) = BV* + APC* , (45) 

BB\ and 



where P solves the Stein equation P = APA* 
the state-covariance S is given by 



(46) 



S = KCPC\ 

where P solves the Lyapunov equation P = APA* + BB*. 

VI. Conclusions and future work 

The ideas and results in [5] were shown to carry over 
to the case where not all interpolation points are at zero. 
This freedom of chosing the interpolation points can be used 
to obtain an improved matching at some frequency regions. 
One example was given to illustrate the effect of moving 
the interpolation points. Input-to-state filters proved to be a 
convenient tool to derive this theory and simple formulas 
based on solving Lyapunov equations were obtained. How- 
ever, if really high order models are considered specialized 
numerical tools have to be developed. 

The approach used in [7] applies only the interpolation on 
S and instead of matching H arbitrary spectral zeros may be 
chosen. This gives the user more freedom in designing the 
model, but at the price of having to tune more parameters. 

In [6], a generalization of the Markov and covariance 
interpolation problem with variable input variance to MIMO 
systems was considered. A similar generalization should be 
possible here. 
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Appendix 

We first derive fijt 

a{I - zA)'^B 



W{z) 



^{I-zA)-^B 
aB + (j{z-^I - A)-^AB 
^B + i{z-^I ~ A)-^AB 
{aB + a{z-^I - Ay^AB) x 

(i^a^-'i-A(3)-'ABm-'] 



B. For state- 


space realizations 











AB{iBY^ ^ 




A state-space reahzation of degree n corresponding to (l38Tl 




-AB£_ 
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AB 




is given by 
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aB J 




W{z) = 


za/3{I - Al3zy^ AB + <7B 


{CB)-\ (43) 




( Ap 
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AB{iB)-^ \ 




where P = I 


— B£^{£^B) ^. From the last section it is known 
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aB j 





that a state-space realization of W of degree n — 1 exists. 




f A/S 


AB 


and one is determined in the appendix. 




\ C7/3 


aB 



Here we have used the matrix 
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to do a change of basis in the large system before cancelHng 
the unreachable second part of the state vector 

This realization is still non-minimal since there is both a 
zero and a pole at infinity. Note that 



I3B ^ B ~ {^B)-^B^B ^ 0. 



(47) 



Now we use F in dTji to do a change of basis. From (|47] | 
it follows that 

cr/3F = al3A [ B ■■■ A"~^B ] , 
T-^A^T = T-^AI3A[ B ■■■ A'^-^B], 



and 



F-^AB = [ 1 •■• 



Then the first state in the new basis is not observable 
or reachable so a reduced order realization is obtained by 
cancelling it: 



W{z) 



where 



and 



A 



TAPAT 


ei 


aPAT 


aB 


' AB ■ 


■ A 





(48) 



and the subindex 2 : n denotes rows 2 to n of the matrix. 

In particular, if the characteristic polynomial XA{t) is 
parameterized as 

XA{t) = t" + Xlt""^ + • ■ • + Xn-lt + Xn, 

the dynamics matrix in (|48] | is 



A 



[ / ] T-^A^AT 

[0 I ]T-^a(t - {^B)-^AB(,f^ 



-Xn~l 

1 -X 

-7 -7«-2 - Xn-l 



where 



X= [ Xn-2 ■■■ Xl ] 

^ r 
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and 7fe = £,A''B/{£,B) for fc = 0, 1, • • ■ , n - 2. 



